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Abstract 

The anisotropic lattice fermion actions are investigated with the one-loop perturbative calcula- 
tions aiming at constructing a formulation for heavy quark with controlled systematic uncertainties. 
For the heavy-light systems at rest the anisotropic lattice with small temporal lattice spacing at 
suppresses the discretization error by a power of atrng for a heavy quark of mass mg. We discuss 
the issue of large discretization errors, which scales as a^mg with a<j the spatial lattice spacing. By 
performing one-loop calculations of the speed-of-light renormalization for several possible lattice 
actions in the limit of at — > 0, we show that one can eliminate the large systematic error on the 
anisotropic lattice. 
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I. INTRODUCTION 



In the heavy quark physics, the lattice simulation of Quantum Chromodynamics (QCD) 
is an indispensable tool to compute hadron masses and matrix elements nonperturbatively 
without introducing model dependence. One of the most important hadron matrix elements 
in the B physics is the B meson decay constant fs, for which a number of lattice calculations 
have been performed so far and the systematic uncertainties are under control at the level 
of 15% accuracy In future, further precise calculation, say better than 5%, is necessary 
to constrain the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements more strictly and to 
search for the signature of new physics. 

One of the dominant uncertainties in the lattice simulation of heavy quark is the system- 
atic error associated with the large heavy quark mass, since the lattice cutoff 1/a available 
with current computer power is not much larger than the heavy quark mass rriQ. A con- 
ventional approach to avoid this problem is to restrict ourselves in the region where the 
systematic error is under control {rriQ ^ 1/a) and to extrapolate to the b quark mass us- 
ing the heavy quark scaling law predicted by the heavy quark effective theory (HQET). 
This is unsatisfactory in order to achieve the 5% accuracy, since the possible error scales as 
(amQ)" {n = 2 for the 0(a)-improved action) and thus grows very quickly toward heavier 
quark masses. The extrapolation to the b quark mass could even amplify the systematic 
uncertainty. 

Another method is the HQET-based approach which includes lattice NRQCD and 
the Fermilab method ^. In this method one considers the lattice action for heavy quark as 
an effective theory valid for large heavy quark masses. The advantage of the HQET-based 
approach is the absence of the large systematic error which scales as {arriQ)^. The price 
one has to pay, on the other hand, is the introduction of a number of terms in the action. 
Their associated coefficients have to be determined by matching the effective theory onto 
the continuum full theory. The matching is usually carried out using perturbation theory, 
which limits the accuracy of the lattice calculation. 

Besides the HQET-based approach, a possible way to control heavy quark discretization 
effects is to consider an anisotropic lattice, where the temporal lattice spacing at is much 
smaller than the spatial one 0,1^ . Since for a heavy meson (or a heavy baryon) at rest the 
large energy scale of order rriQ appears only in the temporal component in the momentum 
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space, one can expect that the systematic error arises as {atrnq)"- and therefore suppressed 
as far as at is small enough. The computational cost is not prohibitive if one keeps the 
spatial lattice spacing relatively large. The problem of the matching of many operators in 
the effective theory does not appear, as the theory is relativistic. 

There is, however, a subtle issue discussed in |2] that for a certain choice of the Wil- 
son term in the spatial direction the systematic error may arise in the combination a^mg 
rather than the expected atmQ and the virtue of the anisotropic lattice is spoiled. With an 
alternative choice the error of order {asmg)"' may be avoided but the unwanted doublers 
become lighter and disturbs the simulation of physical states. The authors of even denied 
the advantage of the anisotropic lattice used for heavy quark based on their observation of 
a^mg-like behavior through radiative corrections. In this paper we discuss this issue fur- 
ther by considering a larger set of 0(a)-improved lattice fermion actions and by performing 
one-loop calculations in the limit — where no a^mQ error remains. 

The appearance of large systematic errors scaling as a^mQ is naively unexpected for the 
following reasons. In the at —>■ limit the only source of the discretization error is the 
spatial derivative in the lattice action. In momentum space, therefore, discretization errors 
scale as a^p with p a typical (spatial) momentum scale in the system, which is of order 
Aqcd for the heavy-light mesons or baryons at rest, and the combination a^mg may not 
appear as the momentum of order mq flows only in the temporal direction. This intuitive 
picture should be correct even after radiative corrections, because the large momentum of 
order mg does not flow into the spatial direction in the momentum space, and therefore the 
discretization error in the spatial lattice derivative cannot accompany the heavy quark mass 
mg. It becomes clearer if one considers the limit l/os <^ mq <^ l/oj, because the spatial 
momentum integral runs up to it /as and thus cannot pick up the larger scale mg. 

Here, .n order to urrderstand the reasorr why the urrexpected a.™«-type error nray appear 
in P, |q|, let us consider the energy- momentum dispersion relation at the tree level. We 
consider the at —>■ limit, and the spatial lattice spacing is also kept small enough such 
that we can neglect the error of O(a^) and higher. For the Wilson- type fermions the inverse 
quark propagator is given as 

mq + ijoPo + i J2 liPi + 7^«« + ^(a^), (1.1) 
where denotes the coefficient in front of the spatial Wilson term as defined in fj2.2|) in the 



next section. The term including the anisotropy ^ = as /at is maintained even in the — 
hmit, because that term could remain when scales as ^. To push up the spatial doubler 
mass in the cutoff scale, r^/^ must be kept finite. 

The energy-momentum dispersion relation becomes 



and thus the error of order a^mg appears unless Tg/^ vanishes, for which the doublers become 
light. Since the term [ts/ Ci^'s^q comes from the cross term of the mass and the Wilson 
terms, the origin of the combination a^mg is nothing to do with the large momentum flow 
into the spatial direction. At the tree level we may consider a set of lattice actions in which 
there is no spatial Wilson term by introducing higher derivative operators to decouple the 
unwanted doublers. This class of actions does not have the problem of the ©(flsmg) error 
at the tree level and may be used for heavy quark even for agrriQ > 1. There are higher 
order terms whose coefficient behaves like agmg, but we neglect them as their contribution 
is O(a^) or higher. 

The problem is, then, whether the nice property of these actions is maintained even with 
radiative corrections. In this paper we perform one-loop perturbative calculation for these 
lattice actions and investigate the mass dependence of the rest mass Mi = E{0) and the 
kinetic mass M2 = {d'^E/dpDpLQ, where -^(p) is the energy of heavy quark on-shell. We 
examine the functional dependence of the speed of light renormalization parameter u, which 
is defined such that the relation Mi = M2 is satisfied. If the one-loop coefficient behaves as 
(ac-mg)", the action suffers from the unwanted heavy quark mass dependent error. Because 
we are interested only in the 0((asmg)") errors, we carry out the one-loop calculation in 
the at —>■ limit, where 0{{atmQ)^) errors vanish. The fermion actions we consider are the 
anisotropic SW (Sheikholeslami and Wohlert) action and some special cases of the 

D234 action p\. We find the latter to be useful for applications to heavy quark systems. 

This paper is organized as follows. In Sectional we define the anisotropic fermion actions 
we consider in this paper, and discuss their tree-level properties. The static limit of those 
actions is considered in Section IIIII The one-loop calculation is then given in Section IIVI 
whose results are presented in Sectional Section EH is devoted to our conclusions. Some 




(1.2) 
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technical details are deferred to the Appendices. 



II. ANISOTROPIC LATTICE FERMION ACTION 

We start with the D234 quark action on the anisotropic lattice ^ given by 

Sd234 = atal'^ij{x)Qij{x) (2.1) 

X 

Q = mo + Y, '^mTmV^ (1 - b^alA^) 

-^«*(E^A«^M + Ec^w^/^^i^M^ +T.^f^d^4K (2-2) 

where ao = at and = {i = 1,2,3) are the temporal and spatial lattice spacings 
respectively, and 

(i^o,i^i) = (1,1^) , (bo^bi) = {bt,bs) , (ro,ri) = (rt,r,) , (2.3) 

(cgw! Cgw) = (cgw ^sw) ' ('^05 di) = {df, ds) . (2.4) 

Note that the lattice spacing in front of the Wilson and clover terms is at, not a^. This 
notation is similar to the one in jsl, but different from those in 0, Q]. The anisotropy 
parameter is defined by 

e = as/at. (2.5) 

The lattice covariant derivatives V^, A^, V^jA^ and A^ represent D^, D^, and D'^, 
respectively, in the continuum theory, and their detailed definitions are given in AppendixIXl 
In this paper we always set 

rt = l, bt = dt = 0. (2.6) 

Thus, the operator Q is nothing but the Wilson-Dirac operator as far as the temporal deriva- 
tives are concerned. With this condition, the energy-momentum relation for the fermion has 
a physical solution only, and the unphysical temporal doublers do not appear jsj. 

Solving Q{p)Q{py = in the momentum space and then setting pq = iE, we obtain the 
energy-momentum relation for the D234 action as 

A ■ u^t^tE. i^^ajEiPi\l + bsa^iPi^)^ + nip f 

4smh — - = — — — , 2.7 

2 l + ywp 
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where 

li{p) = atmo + ^rsttt ^ p^^ + vd^at ^ afpf , (2.8) 
^ i i 

and and p^ are defined in Appendix 1X1 

From (|2.7|) . we obtain the tree-level rest mass Mi = E{0) and kinetic mass M2 = 

{d'E/dpi)-U as 

ajMi = log(l + atmo) , (2.9) 
1 2i^2 

= \ f . (2.10) 

atM2 ajmo(2 + atmo) 1 + ajmo ' 

On the anisotropic lattice with atmo <^ 1, the tree-level mass ratio M1/M2 can be expanded 
in terms of atmo as 

M1/M2 = 1 + (r, - 1) a^mo + 0{{atmof) , (2.11) 

where we set z/ = 1. The deviation of M1/M2 from unity is a lattice discretization error 
arising from the fermion mass. Unless oc C,, such an error is a function of atmo alone, 
wMch is .nail on the anisotropic lattice w.th « 1 rf. When cc J, a d.cteti.ation 
error of order a^mo = ^atmo arises, which is still large on the anisotropic lattice. 

From ()2.7|1 . we can also calculate the "spatial-doubler" mass E'^, i.e. the energy at the 
edge of the Brillouin zone {pi = vr/as): 



atE'' = log 



2rs IGuds 
1 + atmo + -prrid H — Ud 



(2.12) 



= log[l + atm^], (2.13) 

where rid (= 1, 2, 3) is the number of spatial direction with pi = Tr/ag, and the bare spatial- 
doubler mass m^ is given by 

2r 

asmo = a^mo + -^Ud + IGudsUd (2-14) 

in units of the spatial lattice spacing a^. We note that one has to take oc or rf^ > in 
order to decouple the spatial-doubler with the energy E'^ from the physical state for large 
values of ^. 

It is interesting to consider the energy-momentum relation in the Hamiltonian limit at — >■ 
(^ ^ oo), where atmo errors vanish. In this limit the left-hand side of (|2.7|) is replaced by 
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afE , and the energy- momentum relation is simplified to 

E'{a,^0) = ly'Y.P^"i^ + bsa^p")" + ^ (2-15) 

= ml + z/2p2 + (^-^ + 2u%s + 2asmQud)i alY.pt + 0{aY) . (2.16) 

From a small p expansion we obtain a non-relativistic expression of E{at 0): 

E(ai ^ 0) = mo + 7— - ^(p')' + 7^ + 2^^'^^ + Sa^moz/rf, + 0(ay ) . 

(2.17) 

Taking the static limit mg oo subsequently, one obtains 

E{at-^0) mo + iydsalYpt + 0{aY). (2.18) 

i 

Note that the J2iPi term survives even in the static limit. This will be discussed later. 

We are now ready to define our anisotropic actions more explicitly. We study two actions: 
one is the SW action [3,0], and the other is a variant of the D234 action j^, which we call the 
sD34 action. We give these actions and discuss their tree-level properties in the following. 

A. SW action 

The Sheikholeslami-Wohlert (SW) action 0] is defined by 

i^ = r, = 4^ = l, bs = d, = 0. (2.19) 

An 0{at) error arising from the Wilson terms is removed by Cg^ = 1. Since the Wilson 
terms are 0{at), this action goes over to the "naive" quark action in the at ^ oo) 
limit. The energy splitting between the physical state and spatial-doublers E'^ — E vanishes 



in this limit, as one can explicitly find from ()2.14|1 . The energy-momentum relation for the 
SW action is shown in Figure H The energy at the edge of the Brillouin zone decreases as 
^ increases, which shows the reappearance of the spatial-doubler. 

Since = 1, the tree-level mass ratio M1/M2 (Eq. 1)2.111) ) contains no O((asmo)") error: 
M1/M2 = 1 + O((atmo)^) "^^-^ 1. The anisotropic SW action has been applied to the 
simulation of charmonium and the charmed hadrons 11, on ,^ ~ 4 anisotropic 
lattices. 
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B. sD34 action 



We define the sD34 action as 

r, = c^w = , &s > , 4 > 0. (2.20) 

Altliougli the spatial Wilson term is absent (r^ = 0), this action is doubler-free because 
dg > 0. The energy splitting E'^ — E remains even in the at —>■ limit as far as dg is a 
constant independent of ^. Setting (i^ = 1/8 for this action gives the same spatial-doubler 
mass rriQ as for the ^ = 1 SW action when u = 1. The name "sD34" is a reminder that the 
spatial 7i-Df and Df terms survive in the at ^ limit. The sD34 action is similar to the 



one proposed in jl3|, except that those papers consider the case of the isotropic lattice 
^ = 1. Since the sD34 action has the next-nearest neighbor interactions such as ipV^A^ip 
and ipA'jj^tp, this action is more costly to simulate than the SW action which consists of the 
nearest neighbor interactions only. 

The sD34 action does not generate 0{as) discretization errors because of = c^^j = 0. 
In order to remove an 0{at) error arising from the temporal Wilson term with rt = 1, we 
take 

1^ = 1 + ^natmo , = i . (2.21) 

This condition is obtained by performing a field redefinition ipc = ^c'^P, ipc = i^^c, = 
fic = 1 — \rtat{po — rric) to the continuum quark action ipc{x){p + 'mc)i'c{x) 5]. Since 

= the tree-level mass ratio M1/M2 again contains no ©((a^mo)") errors: M1/M2 = 
1 + O{iatmof) 1. 

In the rest of paper we consider the following three choices of bg and dg parameters: 

bs = l, dg = l : sD34 , (2.22) 

o 8 

bg = l, dg = l : sD34(v) , (2.23) 

o o 

bs = l, dg = ^ : sD34(p) . (2.24) 

The first choice sD34, where bg = 1/6, eliminates an O(a^) error arising from the 7jVi 
term: 7iVi(l — ^a^Ai) = 'jiDi + 0(af). The second choice sD34(v), where bg = 1/8, 
eliminates an O(a^) error in the one-gluon vertex ()A13jl : Vi^{q,q\k) = —igt°''~^i + O(a^). 
The difference between the one-loop results for sD34 and for sD34(v) is numerically small as 
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shown in the next section. With the third choice sD34(p), the hopping terms in the action 
are proportional to the projection matrix 1 ± 7^: using the Wilson's projection operator 
^ _ a^7^V^ — |a^A^, the space-component of the action is given by Wi + \wf {i = 1, 2, 3). 



w 



Therefore the third choice can reduce simulation costs compared to the other two choices 



At the tree level the sD34 action ()2.22|1 is 0(ag)-improved, while others contain some 
O(a^) errors. Within the current set of operators ()2.2|) . therefore, the best available choice 
to suppress discretization effects is the sD34 action. 

The energy-momentum relations for the sD34 and sD34(p) actions are shown in Figures El 
and El respectively. For both choices the energy at the edge of the Brillouin zone increases 
as ^ increases, in contrast to the case of the SW action. In small region, the energy- 
momentum relation for the sD34 action is quite close to the continuum one because it has 
no O(agP^) errors. Moreover, in large a^p region near the edge of the Brillouin zone, it is 
close to the continuum one too, for large values of ^. 

To summarize, both the SW action and the sD34 action do not generate the ©((osmg)") 
[n = 1,2, ■ ■ ■) errors at the tree level in the mass ratio M1/M2. While the SW action suffers 
from the spatial doubler for large values of ^, the sD34 action is doubler-free for any value of 
^. Both actions can be used for simulations of the charm quark, if ,^ ^ 2 — 4 and atiric ^ 1. 
But simulations of the bottom quark keeping atrrii, <^ 1 require ^ ~ 5 — 10, for which the 
anisotropic SW action may be contaminated by the spatial doublers. 

Besides the above actions, two other anisotropic actions have been p rop osed and applied 




JJ,UJ, and the 



to heavy quark systems: one is the action with rt = 1 and = ^ p, lli 
other is that with = = ^ [3, ll^ • However, these actions has the spatial Wilson 
term scaling as = ^ and therefore generate the ©((asmo)") errors in the mass ratio M1/M2 
even at the tree level when u = 1, as discussed before and in 7]. For this reason, we do not 
consider these actions further in this paper. 

III. STATIC LIMIT mg ^ 00 AND THE HAMILTONIAN LIMIT at ^ 

In this section we discuss the static limit rriQ — > 00 of the anisotropic fermion actions. At 
finite at, the action always approaches to the usual static action in the limit of atUiQ — > 00. 
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This is shown, e.g. by rescahng the fermion field in ()2.ip as 

„—atMi-t 

i^ix) = ^^hix) , (3.1) 

and then taking ajmo —>■ oo. 

On the other hand, the action in the hmit of mg oo while keeping the condition 
O't^Q ^ 1 can be different. In the at limit, the lattice Dirac operator ()2.2|) becomes 

Q{at ^ 0) = mo + 70^0 + ^T. 7.Vi (1 - 6.a,' A,) + z/^a' E (3-2) 

unless Ts oc ^. Taking subsequently the static limit, the fermion field splits as usual into large 
and small components in the Dirac representation of the Dirac matrix, and the off-diagonal 
terms drop out, the action becomes 

Q{at ^ 0) mo + 70^0 + ^^4a^ ■ (3-3) 

i 

We note that the a^A? term, proportional to vds, remains in the static limit. The SW 
action with dg = ^ approaches the usual static action, but the sD34 action with ds > Q 
does not. This observation is consistent with the static energy evaluated in the — > limit 
fl2.18|) . Formally the static limit ()3.3p can be derived by applying the Foldy-Wouthuysen- 
Tani transformation 



ip{x) — > exp 



E 7iVi (1 - fe.a^Ai) 



i){x) (3.4) 



2mo 

to ()3.2p and then taking mo — >■ 00. 

Results of the one-loop calculation at = in the next section should be consistent with 
the form ()3.3|1 in the mq 00 limit. Suppose that the static action is renormalized as 

(mo + 5m) + 70^0 - \ Srtts E Ai + {uds + 6d)al ^ A^ , (3.5) 

then the mass shift 6m, the kinetic term renormalization 6r, and 6^ do not depend on mo, but 
may depend on udg because the static propagator and vertices contain uds through ()3.3|) . 



IV. ONE-LOOP CALCULATION IN THE HAMILTONIAN LIMIT 

In this section we present the one-loop calculations for the anisotropic actions defined 
in Section IHI We calculate one-loop corrections to the rest mass and the kinetic mass 
renormalization factors in the Hamiltonian limit at 0. From the latter we obtain the 
one-loop correction to the u parameter. 
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A. Formalism 



In the one-loop calculation we basically follow the notation of 21| and add some exten- 
sions to the case of the anisotropic lattice. 

We write the inverse free quark propagator as 

atG^\p) = iatI^{p) + a,L{p), (4.1) 

and the self energy as 

atS(p) = i 7m^m(p) sin(a^p^) + C{p) 

^ iY.^,B,{p) + C{p), (4.2) 

where p is the external momentum. The inverse full quark propagator is then given by 

G-\p) = G,\p)-np)- (4.3) 
Solving G^^{G^^y = with pq = iE, we obtain the all-orders dispersion relation 



1 + /x(p) - cosh(atE) -G=J{1- smh\atE) - EAatKi - sm{a,Pi)Y , (4.4) 



where /i(p) is given in 

Setting p = in (|4.4p . we obtain the rest mass Mi = E{0) as 

gatMi ^ l + atmo + Ao{iMi,0)smh{atMi) -G{iMi,0) 

= l + atmo-iBo{iMi,0) -G{iMi,0). (4.5) 

In order to have massless quarks remain massless at the quantum level, we need a mass 
subtraction. Defining the critical bare mass atrriQc = C(0, 0), we can write 



1 + atMo - iBSMi, 0) - aub(^Mi, 0), (4.6) 



where Mq = mo — and Csub(^^i, 0) = G{iMi, 0) — atniQc- When Mq = 0, the rest mass 
Ml vanishes by construction. Usually the mass subtraction is done nonperturbatively in the 
numerical simulation by defining the critical hopping parameter. 
In perturbation theory the rest mass is expanded as 



oo 



Mi=5:/'Mi^ (4.7) 
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in which the tree-level rest mass is 



asM? = ^log(l + atMo) "-^ a,Mo, (4.8) 

while the one-loop coefficient is given by 

asMP = (-^^B^^\^M^, 0) - ^C^d^M,, 0)) e-"'<' 

-^^B^\^Mu 0) - ^C^Sdm, 0). (4.9) 

Note that Mi is now normalized by the spatial lattice spacing a^. Before the subtraction 
the one-loop coefficient is given by 

«Xlosub -^ei?fl(^Mi,0)-ecW(^Mi,0). (4.10) 

Differentiating ()4.4p in terms of pi twice and then setting p = 0, we obtain the kinetic 

mass M2 = {d'^E/dpDplg as 

eM _ A,(iM,, 0) cosh(a,MO _ r, [u/^ - A^jtM,, 0)]^ 

eatM2 e ^ ^ [1 - Ao(2Mi, 0)] sinh(aiMi) ' ^ ' ' 

where 

^(0) = ;77^^^[^o(^i?(p),p)sinh(a,Mi)-C(zE(p),p)]p^o 



with 

DUO) 



92 ri 

-Bo-C 



d{asPiY 
d_ 

d{asPo) 



p=(iMi,0) 



(4.13) 



The kinetic mass is expanded as 



M2 = Y.9''m!\ (4.15) 

1=0 



and the tree level relation becomes 



0*^2 sinh(aiMj 

where iVff = m2(M{°'). 



roK (4.16) 
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From ()4.11|) and ()4.16|) . we can obtain the kinetic mass renormalization factor defined by 



M. 



'M2 



' -l + E/'<- (4-17) 



"^2 (Ml) 

Here the argument of m2 is the all-orders rest mass Mi. The one-loop coefficient is given by 

[1] ^ 2piA^l\iMu 0) - z/MgJ(^Mi, 0) - DM(0)e^ sinh(a,Mi) 

z/2 + r, sinh(atMi) 

-A!)^'(2Mi, 0) cosh(a4Mi)e-"*^^\ (4.18) 
In the at — hmit, the tree-level kinetic mass goes to 

— '^°^ + i?.a,, (4.19) 

m2 Ml 



where we defined Rg = r^/^, and hence goes to 

Zm, + i?.a,M2 . (4.20) 

Therefore, in this limit, the renormalized u parameter and Rg which give Mi = M2 can be 
determined from Zm2'- 

= 1 + (2z/W + RfasMf^)g'' + 0{g^) (4.21) 

with the one-loop coefficient 

<=2z.W+i?Wa3Mf . (4.22) 

Here we used "^^^ 0. 

On the other hand, from ()4.18p we obtain in this limit 

Z'l - e4^^(^Mi, 0) - ^ ^"^"'^^^^ °^ - a.m,ei^W(0) , (4.23) 



where A^l\iMi,Q) = |^|p=(iAfi,o). 



9(aspi) I 

B. One-loop diagrams 

Here we compute one-loop contributions relevant to the rest mass and the kinetic mass 
renormalization. At the one-loop level, the self-energy is written as 

S[i](p) = S^^^(p) + S*^^(p) + ST-i-(p), (4.24) 

13 



where the contribution from the regular graph Fig. EJ^a) is denoted by H^'^^lp), while the 
tadpole graph Fig. m^b) gives S*^*^. In order to remove the bulk of S****^, we apply the tadpole 



22], which amounts to E"""'^'. Feynman rules relevant to our calculations are 



improvement 

summarized in Appendix^ We use the anisotropic Wilson gluon action given by 



S -1 
'~ 9' 



x,i>j 



(4.25) 



where Poi{x) and Pij{x) are the temporal and spatial plaquettes respectively. In the calcu- 
lations of Y7°^ and T}^, we adapt the Feynman gauge a = 1 for the gluon propagator. 

C. regular graph 

The contribution from the regular graph is 



1 



. (27r)4 ^2^2 ^ ^2 K^ip - kf + L{p - k) 



2 ' 



(4.26) 



where L = atL, K = atK and the gluon momenta are rescaled as a^k^ — ^ k^ and a^k^ 
and 



p 

Fc = atFc = 2Y.K,X,Y,-LY.{Xl-Y^ 



- p ) 



(4.27) 
(4.28) 



with K = K{p -k), L = L{p - k), = X^{2p - k, ±k) and = Y^{2p - k, ±k) given in 
Appendix 1X1 Since the vertex from the clover term (Jij,yF^u in the fermion actions are 0{at) 
and vanishing in the — > limit, we omit their contributions. 
In the — i> limit, where 



'BO 



(Po-A:o){l-E(^| + >^')}> 

i 

j 

2J:K,X,Y,-L{1 + Y.{X]-Y^')}, 



(4.29) 
(4.30) 
(4.31) 



14 



we obtain 



e^rb) [ 2,2\ 12 k) S,{p - k) , (4.32) 

eC-s(p) / ^ ^tAt^ Fc{v, k) S,{p - A:) . (4.33) 



CO 



and 



Here we defined 

^ (Po-*o)^ + aAi(p-*)^ + i(p-*)^- 

Differentiating ()4.32p and ()4.33p in terms of external momenta p and then setting p = 
(iMi,0), we obtain tlie contributions to m|^^ and Z^^}^ according to ()4.9p and ()4.23p . For 
tlie evaluation of the loop integrals, we first integrate over /cq analytically as described in 
Appendix ^ The remaining integration over k is evaluated numerically using an adaptive 
integration routine VEGAS [2^. 

Since the rest mass and the kinetic mass are physical quantities, one-loop corrections to 
them are infrared-finite. Although there are infrared divergences in the partial derivatives 
Dis{0) and Dit{0) in the kinetic mass renormalization, they cancel in the total derivative 
D{0). In numerical integrations over k, we evaluate the total derivative directly, rather than 
evaluate each partial derivative with subtraction of the infrared divergences. 



D. tadpole graph 

Although the calculation of S****^ is much simpler than that of E^''^^, it is worthwhile to 
show the dependence of the results on the mass and on the parameters. The contribution 
from the tadpole graph at finite at is given by 



as^^'^'ip) 



. (271)4^2^2 + 



at 



i'y^ I a^X^ sin(a^p^) + 4a^Z^ sin(2a^p^) cos^(^) 



Xfi cos(a^p^) + Aa^W^ cos(2a^p^J cos 



k 



(4.36) 
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from which we immediately obtain 



^B';^ip) = Cf^ {Ja^X^ + 8T^ a^Z^ cos(a^p^)} sin(a^p^), (4.37) 
iC'^'ip) = -CfY. — {-^^m^m cos(a^p^) + 4 a^W^ cos(2a^p^)} , (4.38) 



where 



J = U ^ . ^ . , (4.39) 
^ (27r)4^2^2 + k2' ^ ^ 



T 



^1-. (27r)4 ^2^2 ^ k2 "^^^ ^y^- ^- 
In the tti hmit, J = Tq = 0.2277 and = 0.1282. 

Tadpole contributions to m|^^ and Z^^j^ in the at —>■ limit are easily calculated from 
()4.37|1 and ()4.38j) . The contribution to the rest mass before the subtraction ()4.10|) is given 

by 

«.M^;^nosub C^^{Ja,F, + 4T,a,lV,}, (4.41) 

i 

which depends on agYi and agWi , z.e. udg, but not on the mass. It contributes to the critical 
mass a^moc only, so a.sMl^'^ = after the subtraction. 

The contribution to the kinetic mass renormalization is given by 



2 

- {JagXi + 8 Ti asZi) + asm2( Ja^Fi + 16 Ti OsW^i^ 



(4.42) 



We find that a term proportional to as"^2 appears, which depends on udg again. This 
manifest asmq dependence originates from the ^-D(is) term in ()4.23|) . Therefore diverges 
as Olasniq) toward the static limit for the sD34 action with dg > 0, while it is mass- 
independent for the SW action with d^ = 0. Similar mass-dependences are also observed in 
^mI- We will discuss the ©(a^mg) divergence of Zm2 in Section IVl 



E. tadpole improvement 

Tadpole improvement 



22j | is achieved by replacing the link valuable by U^/u^^ where 
Ufj_ = (U^) is the mean link valuable. In perturbation theory the contribution from the 
tadpole improvement is obtained from the difference between the inverse free propagator 
Gq^ and the tadpole-improved inverse free propagator (G^^)'^'^' |18ii]. In momentum space 
the latter is given by the former with replacements 

sm{na^p^) ^ sm{na^p^)/u'^, cos{naf,Pf,) cos{na^p^)/u'^, (4.43) 
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where n = 1, 2, • • •. The one-loop contribution is then given by 

a.ST-i-(p) = {a,G,\p)-a,iG,Y-'ip))/g' 

= [^^M {2a^X^ + 8a^Z^ cos(a^p^)} sin(a^p^) 

- {2a^F^ cos(a^p^) + 4a^W^ cos(2a^p^)}] , (4.44) 

where we expanded 

«,= l + /nW+0(^7^) (nW<0). (4.45) 
We adopt the mean hnk in Landau gauge for the definition of u^, which is given by 

«W = -l^C, j;=° (4.46) 

^ CI Q 



with 



. , . , . ^ , . . , (4 47) 

- (27r)4 ^2^2 + k2 1 C^o+k^ ' 

= J - (1 - a) 5J^. (4.48) 



In the — > hmit, 5Ji = 6Jo = \ J ■ We then obtain 



iBl^iv) = -Cf^J^=' {a^X^ + Aa^Z^cosia^p^)}sm{a^p^) , (4.49) 

Of. 

^ as 

Comparing and (p3n|l with (jOTjl and (jOHI), we find that and ^Cj-^- are also 



^C^-'ip) = CfY. —Jr' cos(a^p^) + 2a^W^ cos(2a^p^)} . (4.50) 



obtained from ^-Bjf^ and ^C*^'^ with replacements 



^ - - ' ^>^^-l ^r°- (4-51) 

Contributions to Mj^j^Qg^^^ and Z^l^^ from the tadpole improvement are given by Eqs. ()4.41|l 
and ()4.42p with the above replacements. Since ~ f"^' tadpole contributions are largely 
canceled by the tadpole improvement. 



V. ONE-LOOP RESULTS 
A. Rest mass 

Now we present the results of our one-loop calculations in the at ^ limit. The one- loop 
correction to the rest mass a^Mj^' is plotted as a function of agMf'^ / {1 + agMi'^) in Figure|3l 
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and numerical values of OsM]^' and asmo). are given in TableHl As shown in the figure, 
for all the actions increase from the massless limit and reach their maximum values around 
agMi'^ = 1-3, then decrease to the static values represented by open symbols. Fitting our 
results in the small mass region (a^ mP < 1), we confirmed that the one-loop corrections 
are consistent with the mass singularity 



Mf 1 ~ -C^FT^^Mf 1 log(a.Mf 1)2. (5.1) 



The results for a^Afj^l in the static limit agMf'^ oo depend on the action, since the 
reduced static action ()3.3|1 includes the a^Af term proportional to udg. This situation is 
in contrast to the case of finite at calculations in 18[ 21 1, where cts^i^^ goes to a universal 



value. In the finite at case, a^Mj^l does not depend on udg in the a^mg — > oo limit, because 
the static action always gives the Wilson line. 

B. Kinetic mass renormalization 

The one-loop correction to the kinetic mass renormalization is related to the speed 
of light renormalization u according to ()4.22j) . The study of the agniQ dependence of u at 
the one-loop level is a main purpose of this paper. The results of Z^^j^ are shown in Figure El 
and their numerical values are given in Table lUl 

First, we focus on the result for the SW action, which becomes the naive quark action 
in the — > limit as the Wilson term and the clover term vanish. From Figure El (lower 
panel), we find that the mass dependence of Z^^j^ (filled circle) is very weak, and Z^^}^ stays 
constant in the infinite mass limit. A difference between the value in the static limit and 



0) = —0.006. This is only 6% of the same 
2l|. The result implies that mass dependent 



that in the massless limit is Z^l}^{oo) — Z^^^ 
difference for the isotropic SW action —0.10 
discretization errors of order g^iagrriQ)"' for are small on the anisotropic lattice. The 
same conclusion holds for any action which becomes the naive quark action in the at ^ 
limit. For instance, the action with = and dg = d/C,, where is a constant independent 
of C,, belongs to this class. However, we remark that such actions suffer from the spatial 
doublers for large values of ^, as mentioned in Section ITTl 

Next, we consider the results for the sD34 actions, which are doubler-free even in the 
at ^ limit. As shown in Figure IHl (lower panel), Z^^j^ for the sD34 actions monotonically 
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decreases as the mass increases, and diverges as Oi^asmq) toward the static hmit. 

The 0{asmQ) divergence of Z^l}^ is due to the finiteness of 'C-D(is) in the static hmit 
multiphed by asm2 in (I4.23|) . The appearance of this manifest asfng dependence, which is 
proportional to ud^, can be explained as follows. Using Z]^}^, the kinetic term renormalization 
6r for the static action (|3.5|) is given by 



1 1 \ . o Z 



[1] 



6r= lim (_^-^)=- hm g^^ + Oig'). (5.2) 

rriQ^oo \ asM2 asm2/ mg^oo asm2 

Because 6r is a constant independent of the mass, zj^^ diverges as 0{asmQ) in the large 
mass limit. 

However, we note that this kind of O {asfng) divergence is nothing to do with the dis- 
cretization error increasing as agfriQ but a renormalization of the reduced static action ()3.5p . 

"A/2 

ziLn^ = zil + 6l}^ a^m^ (5.3) 



In order to isolate such an 0{g^asmQ) effect, we consider a subtracted Z^^}^ defined through 



as a measure of the remaining 0{g'^{asmQ)^) {n > 2) errors. After the subtraction of the 
manifest a^mQ dependence, Z^^j^ for the sD34 actions converges to a finite value in the 
static limit as shown in Figure [7| We also find that the mass dependence of Z^^^ sub f*^^ 
sD34 and sD34(v) actions is as small as that for the SW action. Note that Z^^j^ = Z^^^ 
for the SW action because of 6r = 0. 

Another way to discuss the remaining 0{g'^{asmQ)^) {n > 2) errors for the sD34 actions 
is to assess their linearity in the mass parameter. Since z]^}^ for the sD34 actions seems 
like a linear function of a^Ml^' effectively as shown in Figure El (lower panel), we attempt a 
linear fit using the data for a.M^ < 0.5. The fitting hues = Z^^l{0) + x a.iVff 

shown by dashed or dotted lines approximate Z^lj^ very well from the small mass region 
a.sMf'^ ^ 1 to a relatively large mass regime a^Mj^' ~ 1. The difference Z^^j^ — Z^lj^^m 
plotted in the upper panel of Figure El We find that the difference is less than or about 
0.005 (0.01) at agM^^ = 1 (3) for the sD34 and sD34(v) actions, and slightly larger for the 
sD34(p) action. Since the (renormalized) coupling constant is g^ = Anas ~ 2 in current 
simulations, the difference from the linearity g'^{z]lj^ — Z^lj^ ^^^) is small compared to the 
tree-level value zfj^ = 1. It indicates that 0{g'^{asmQ)^) {n > 2) errors are suppressed on 
the anisotropic lattice, and the sD34 actions can be well approximated by a linear 

ansatz; z[}^ ^ ^i/a.iin- 
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If one would like to avoid the appearance of the renormalization scaling as agmq, it is 
possible to tune the spatial Wilson term as = — (jW such that the second term in ()5.3p 
vanishes, and then the one-loop coefficient of the speed-of-light renormalization is given by 
z/t-*^! = z]^j^ ^^■^^/2. Since the remaining ©((asmg)") correction for u^'^^ is small and does not 
diverge as a function of agrriQ as shown in Figure [7[ it essentially solves the problem of large 
radiative correction in the anisotropic lattice actions for heavy quark. It also suggests that 
if one can nonperturbatively tune the Wilson term in the static limit, e.g. by adjusting 
until the ©(a^mg) divergence of for mesons goes away, the above cancellation of the 
asTTiQ error can be implemented nonperturbatively. 

VI. CONCLUSIONS 

In this paper we discuss on the issue whether the discretization error scales as (asmg)" 
when the heavy quark action is discretized on an anisotropic lattice for which the temporal 
lattice spacing at is very small in order to keep the condition atrriQ ^ 1 while the spatial 
lattice spacing is relatively large and agmq can be order one. Our naive expectation is 
that the discretization error does not behave as a^mg for the heavy-light mesons (or baryons) 
at rest, since momentum scale flowing into the spatial direction is of order of the QCD scale 
Aqcd rather than the heavy quark mass scale mg. Even at the quantum level the maximum 
(virtual) momentum flowing into the spatial direction is vr/a^, and the discretization error 
coming from the spatial derivative cannot pick up the large heavy quark mass. 

Through the one-loop calculations of the kinetic mass renormalization for a class of 
lattice fermion actions, we found that our expectation is indeed the case. For the sD34 
actions there is a piece which behaves as agmq in the one-loop coefficient of the kinetic 
mass renormalization, but it originates from the renormalization of the spatial Wilson term, 
which remains even in the static limit, and thus does not come from the discretization of the 
spatial derivative. It implies that if one can nonperturbatively tune the spatial Wilson term 
(the parameter r^) such that it vanishes in the static limit, the unwanted behavior a^mq can 
be removed from the speed-of-light renormalization. Although there is a possibility that the 
unwanted discretization error scaling as agmq exists in some other quantities, it is unlikely 
from our considerations. 

The anisotropic lattice thus remains as a promising approach to treat heavy quarks on the 
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lattice. As in the usual relativistic approach, the theory is renormalizable and the number 
of necessary terms in the action is limited. It also opens a possibility to tune the parameters 
in the action nonperturbatively for heavy quarks. 
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APPENDIX A: DEFINITIONS AND FEYNMAN RULES 

The lattice covariant derivatives are defined by 



1 

2a 



1 



U^{x)i!{x + fj.) - U_^{x)ij{x - fi) , (Al) 
U^,{x)ij{x + fj,) + U^^{x)7jj{x- fj,)-2^{x) , (A2) 
U^{x)U^{x + ij)ii{x + 2//) - C/_^(x)C/_^(a; - iJi)ii{x - 2/1) 



-2U^{x)^{x + 11) + 2U_^{x)jP{x - ii) 



U^ix)U^{x + fx)ij{x + 2/x) + U_f,{x)U_^{x - fx)^{x - 2/x) 



AU^{x)ij{x + n)- AU^^{x)i^{x -n) + 6V^( 



x\ 



(A3) 



(A4) 



We also define the lattice momenta 



a^p^ = 2sin(a^p^/2). 



(A5) 
(A6) 



Feynman rules for our anisotropic actions can be derived in usual way. The gluon prop- 
agator with Feynman gauge is given by 



A;2 



(A7) 
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The quark propagator is 



where 



Koip) = Po Po, (A9) 

K,{p) = up,{l + hsalPl), (AlO) 



and 



L{p) = mo + ^atY,r^pl + ud,Y,a^p^ (All) 
m, + udsY.4pt (A12) 

i 

for our quark actions with Eq. ()2.6|) . 

The one-gluon vertex with the incoming quark momentum g, the outgoing quark momen- 
tum q' and the incoming gluon momentum k = q' — q is given by 

Vi^,^{q, q', k) = -tgt^ [7^X^(5 + q', k) - iY,{q + g', A;)] , (A13) 

where 

X^,{q + g', k) = 2a^X^ cos ^ ^a^^m + Qa^^m ^ _^ 4a^Z^ cos (a^q^ + a^g|,) cos (A14) 
Y,{q + g', fc) = 2a,Y, sin [^^^l^fi±^^ + ^a,W, sin (a,g, + a,g;) cos {^-^^ ( A15) 

and 

a^X^ = 1 1/^ + 1/^5^, (A16) 

a^y;. = h4z/^d^, (A17) 

a^Z^ = -^u^b^, (A18) 
a^,W^ = -i^i^d^. (A19) 

The t"" are generators of color SU(3). We ignore the one-gluon vertex arising from the clover 
terms because such a vertex becomes irrelevant in the at —>■ limit. 

Finally the two-gluon vertex with the incoming gluon momenta k and k' {k + k' = q' — q) 
is given by 

V^;i^{q,q\k,k')=2a,g\tH')x 
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<^ a^X^ sm I I + Aa^Z^ sm [a^q^ + a^g^j cos I I cos I I \ 

- \a,Y,cos [^J!St±^\^ + 4a,W^,cos {a,q, + a^g^) cos cos |(|\20) 

Here we omit terms that vanish by symmetrizing between two gluons and that arise from 
the clover terms, which are unnecessary in the calculation of the tadpole graph. 



APPENDIX B: /sq-INTEGRATIONS 

In this Appendix we summarize some formula on the /co-integrations, which are needed 
for the calculation of the regular graph. We use the following results for one-dimensional 
integrations: 



oo ]^ 

dx 



oo a + hx"^ g'^{ie — x)^ + 



.1 1 11. 

fd) -g^{e- + P a-h{e + f/gf fg^ ' 



f°° , 1 die — x) 

/ dx ■ 

J —ex. 



a + hx"^ g'^{ie — x)"^ + p 

.1 ^~ 1 1 

^ ^ -g2^e - v^)2 + p~ a-h{e + f / gf 7^ ' ^^^^ 

/oo 1 c 1 dli 

77¥77— = -TT71^' (^3) 



h 
h 



7 



1 


c 


1 dh 


a + bx'^ {g'^{ie 


-a;)2 + /2)2 


2/ 9/' 


1 c( 


ie — x)^ 


1 dh 


a + bx'^ {g'^{ie 


-a;)2 + /2)2 


2g dg' 


1 c 


— x) 


1 dh 


a + bx"^ {g'^{ie 


-x)2 + /2)2 


2/ 9/' 


1 


C 


1 9/3 


a + bx"^ {g'^{ie 


-a;)2 + /2)3 


4/ 9/' 


1 c 


ie — x) 


1 a/5 



= rdx^ dte-xr ^ _^dh . . 

~ i-oo a + 6x2 ( 2(i^ _x)2^f2)2 2 f 9 f ^ ^ 



r , _^ C(2e - X) ^ _ J_dh (r,n. 

y_oo a + 6x2 (^2(^g _ ^)2 + j2)3 4j ' y ) 

where e < f/g is assumed. These integrations are calculated by hand using the residue 
theorem, and checked by Mathematica. 

In the calculation of the regular graph, we assign 

x^ko, a^|k^ b^al, g^l, e ^ M, , f ^ E{k) , (B8) 
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where 



E(k) ^ liy^^kf(l + bsa^kfr + (mo + udsj:alkfy 



(B9) 



The overall factors c depend on the spatial momentum k. Using integrations /1-/7 with 
above assignments, relevant contributions from the regular graph are given by 



e^r/(o) 



k ZTT 

k ZTT 



Al + h-Al) , 



1 

k 2^ 



-{h-Ds + h-Ds + h-Ds) — {h-Ds + h-Ds + 



— / —{h-Dt + h-Dt — ih-Dt) ■ 

Jk ZTT 



(BIO) 

(Bll) 
(B12) 
(B13) 
(B14) 
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SW 


sD34 


sD34(v) 


sD34(p) 


0.0 


OOOOOOfOO) 


000000 fOO) 


OOOOOOfOO) 


000000 f 00) 


0.1 


0.033586(11) 


0.018299(19) 


0.018360(21) 


0.015615(23) 


0.2 


0.059382(20) 


0.029235(21) 


0.029306(21) 


0.024206(23) 


0.3 


0.081252(22) 


0.037292(24) 


0.037354(23) 


0.030177(41) 


0.4 


0.100023(23) 


0.043466(28) 


0.043405(21) 


0.034594(25) 


0.5 


0.116008(27) 


0.048228(29) 


0.048196(25) 


0.037909(25) 


0.6 


129786fl9) 


051966f24) 


051762f25) 


040428 f 32) 


0.7 


0.141467(36) 


0.054948(34) 


0.054572(31) 


0.042396(31) 


0.8 


0.151341(22) 


0.057183(24) 


0.056763(22) 


0.043836(33) 


0.9 


0.159679(29) 


0.059029(36) 


0.058465(29) 


0.044916(49) 


1.0 


0.166758(21) 


0.060338(33) 


0.059687(24) 


0.045724(34) 


2.0 


0.196227(24) 


0.062578(35) 


0.060997(28) 


0.045709(32) 


3.0 


199166('39') 


058705(39) 


056713('28) 


041241 f34) 


4.0 


0.197161(26) 


0.054263(23) 


0.052280(29) 


0.036472(30) 


5.0 


0.194441(44) 


0.050519(28) 


0.048617(55) 


0.032255(53) 


10.0 


0.184712(29) 


0.038993(46) 


0.037121(28) 


0.018373(37) 


oo 


0.168490(26) 


0.019045(21) 


0.018249(29) 


-0.009803(27) 


[1] 




0.000000(00) 


-0.060828(02) 


-0.061672(02) 


0.001431(03) 



TABLE I: Numerical values of a^Mj^ for various values of a^Mj , and asmpj for the SW action 
and the sD34 actions. 
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n n/i 7QOfi c;n^ 
— u.U4/yyDi^ouj 


U. 1 


U.Ui ( OU / \ J 


— U.UUooo t yoZ j 


n nni qr^/'^9~i 
— u.uuiyoo^^oz J 




U.O 


U.Ui ( OZO(^0^ ) 


n ni 99QQi'S/i^ 
— u.uiz/yy (^o4 J 


— U.UU004o(^Zo j 


n nfii^^nni' ^n^i 

— U.UDOOUU^^OU ) 


n Q 
u.y 


n ni 71 Aot A7\ 

U.U± 1 ±4Z^^4i ( J 


n ni R9'^S('9s'i 

— U.UlDZoo(^Zo J 


n nnoi 1 7/'9s^ 

— U.UUyil / yZo ) 


n n7/in7n('/iR'i 

— u.u ( 4U / U(^41:D ) 


l.U 


n ni fiSfio/'^fi'i 

U.UIDOOZ^^OD J 


— u.uzui^y \^zo j 


n ni 9fifii \ 

— U.UlZODH^Ol J 


— U.U(5Z40D\^00 j 


z.u 


n ni 1^1 7f{('it^\ 


n nl^f^7/l9^'^/l^ 

— U.UOD 1 tZyot ) 


— U.U4iDDD0^^Zo J 


— u.iDuzoy^^oy j 


3.0 


0.014246(21) 


-0.090607(29) 


-0.078766(30) 


-0.229675(78) 


4.0 


0.013677(17) 


-0.122890(42) 


-0.109982(40) 


-0.293550(90) 


5.0 


0.013371(18) 


-0.154177(45) 


-0.140336(48) 


-0.35409(12) 


10.0 


0.012716(15) 


-0.304675(91) 


-0.288888(86) 


-0.63276(22) 


oo 


0.012316(07) 


— oo 


— oo 


— oo 



TABLE II: Numerical values of ZU for the SW action and the sD34 actions. 
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FIG. 1: Energy-momentum relation at different values of ^ for the SW action. The left panel 
shows the case of asE{0) = 0, while the right shows asE{0) = 1. The spatial momentum p is 
along the (1,1,0) direction. For comparison we also plot the energy-momentum relation in the 
continuum. 




FIG. 2: Energy-momentum relation for the sD34 action. 
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FIG. 3: Energy-momentum relation for the sD34(p) action. 




(a) (b) 



FIG. 4: Feynman graphs relevant for the one-loop quark self energy. The left (a) is the regular 
graph, and the right (b) is the tadpole graph. 
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FIG. 5: OsMf' versus agUf'^ /{I + agMf'^) for the SW action and the sD34 actions. The values in 
the static hmit are denoted by open symbols. 
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FIG. 6: The lower figure shows versus agM^^ for the SW action and the sD34 actions. 

The value in the static limit for the SW action is denoted by open circle. Lines are the linear 
approximations to the results for the sD34 actions (-^^^ jin) as explained in the text. The upper 
figure shows the difference — versus agMf^ for the sD34 actions. 
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FIG. 7: ^M2,sub for the sD34 actions together with .^M2,sub = -^Ma foi" ^^^e SW action. 
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